Chaos and Thermalization in a Dynamical Model of Two Interacting Particles 
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A quantum dynamical model of two interacting spins, with chaotic and regular components, 
is investigated using a finite two-particles symmetrized basis. Chaotic eigenstates give rise to an 
equilibrium occupation number distribution in close agreement with the Bose-Einstein distribution 
despite the small number of particles (n — 2). However, the corresponding temperature differs from 
that derived from the standard Canonical Ensemble. On the other side, an acceptable agreement 
with the latter is restored by artificially randomizing the model. Different definitions of temperature 
are then discussed and compared . 

PACS numbers: 05.45. +b 
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The concept of quantum ergodicity, and its connections 
with the foundations of statistical mechanics, have a long 
history which, in the last years, has received new impulse 
from the investigation of quantum systems which are 
chaotic in the classical limit. The important issue here is 
twofold. First, one would like to know if such quantum 
systems display some sort of chaotic behaviour; second, 
if this behaviour will provide foundations for quantum 
statistical mechanics. In other words, the problem is 
whether a N- body isolated system, which, on the clas- 
sical side, displays sufficiently strong chaotic properties, 
will, on the quantum side, give rise to conventional sta- 
tistical properties, such as, e.g., the Bose-Einstein (BE) 
or Fermi-Dirac (FD) distributions 0. The first problem 
has been widely investigated for single-particle systems; 
not so for many-body systems, essentially because such 
systems are not easily accessible to numerical investiga- 
tions. Instead, little is known about the second problem. 

Whereas early studies on the foundations of quantum 
statistical mechanics did not attach special importance to 
the structure of eigenstates, quantum ergodicity is nowa- 
days associated with eigenstates being homogeneously 
spread, in an appropriate statistical sense, over the whole 
energetically allowed range. Much of our current un- 
derstanding of the structure of eigenstates for classically 
chaotic systems is based on the study of single-particle 
systems, and on the analysis of eigenfunctions of Ran- 
dom Matrices, which have been often used as models 
for quantum chaos. In particular, suitably constructed 
Wigner Banded Random Matrices (WBRM) g] are con- 
jectured to preserve many of the spectral properties of 
"chaotic" Hamiltonian matrices H. For WBRM, several 
results are known [^-Q, concerning spectrum statistics, 
structure of eigenstates and of the Local Spectral Density 
of States (also called Strength Function), and conditions 



for quantum ergodicity have been given j7j . These results 
provide paradigms of quantum ergodic behaviour, which 
should be tested on realistic Hamiltonians. 

If quantum ergodicity can be established along such 
lines for quantum many-body systems, then the problem 
arises, whether it gives rise to some sort of statistical 
equilibrium. One would like to know whether the quan- 
tum averages of the occupation numbers of single-particle 
energy levels, taken over many-body ergodic eigenstates 
(or, more properly, over mixtures of eigenstates with en- 
ergies lying in a narrow selected range), yield some sta- 
tistically stable distribution law; if this is the case, how 
large must the number of particles be, in order that this 
distribution reproduces the conventional Fermi-Dirac or 
Bose-Einstein statistics? 

This theoretical approach has been advocated in Refs. 
|p|-fi0[; in this Letter we numerically investigate the just 
sketched theoretical issues on an isolated system of two 
interacting particles, which, under appropriate condi- 
tions, is classically chaotic in some energy region. This 
model has been proposed and investigated in within 
the framework of quantum chaos; here we analyze it from 
the standpoint of statistical mechanics, and show that 
thermalization occurs in the classically chaotic energy re- 
gion, in a sense that will be discussed below. 

We first review some fundamental facts about its clas- 
sical and quantum behavior. The model describes two 
coupled rotators, with angular momentum L , M and 
Hamiltonian : 



H = H + V={L Z + M z ) + L X M X 



(1) 



It may be used to describe the interaction of quasi- 
spins in nuclear physics. Constants of motion are H = E, 
L 2 and M 2 . It is worth to mention that in this form the 
dynamical variables L, M are not canonical. 
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The analysys of the surfaces of section reveals a large 
number of regular trajectories covering invariant tori 
when L 2 ,M 2 are both very small or very large p[ . To 
simplify the problem we set L = M. In such a case 
the most interesting situation occurs when 1 < L < 10 
where, depending on the energy value E, regular and 
chaotic regions coexist. Typically when \E\ is close to 
the maximum allowed energy E max = L 2 + 1 JTl| trajec- 
tories are regular while for E ~ the islands of stability 
become very small and chaotic motion dominates. 

Quantization follows standard rules, and angular mo- 
menta are quantized according to the relations L 2 — 
M 2 = h 2 l(l + l) where I is an integer number. Therefore, 
for given I the Hamiltonian is a finite matrix, and the 
semiclassical limit is recovered in the limit I — > oo and 
h — > keeping L 2 constant. 

The matrix elements in the basis \l z ,m z ) have the 
form, 

(l' z ,m' z \H \l z ,m z ) = S m ^ m 'Ji z ^h(l z + m z ) 

h 2 

(l' z ,m' z \V\l z ,m z ) = -^8 m ^ m > z±1 8i^ Vz± i x (2) 
x[(Z + l z ){l -l z + l)(l + m z )(l -m z + I)] 1 ' 2 

with l z ,m z integers, —I < l z ,m z < I . 

The z-component of the total angular momentum 
J z = L z + M z (which is the unperturbed Hamiltonian 
Hq) obeys the selection rules A J z = 0, ±2h, so the sub- 
space spanned by the states with odd J z can be separated 
from that with J z even (there are no matrix elements for 
the transition between them). In what follows, we fix 
J z = Hq even (multiple of h). 

A key point in our approach is to represent the Hamil- 
tonian in the symmetrized two-particle basis of non- 
interacting particles. This corresponds to the well-known 
"shell model" representation used in atomic and nuclear 
physics. Here, we restrict our considerations only to sym- 
metric states with respect to the exchange of the two 
particles. In the symmetrized basis each set of states 
with fixed even Ho has a degeneracy l + l— \Hq\/2K, and 
the dimension of the Hamiltonian matrix is TV = (/ + l) 2 . 
Finally, we reorder the matrix according to increasing un- 
perturbed energy, and thus obtain a band matrix, with 
high sparsity within the band (each line has at most 5 
elements). 

Direct diagonalization of such a matrix gives the eigen- 
functions tp n (E m ) of the total Hamiltonian H represented 
in the ordered symmetrized two-particle basis |n). Here 
ipniEm) is the n-th component of the eigenfunction hav- 
ing E m as eigenvalue. A detailed analysis of the struc- 
ture of eigenstates, (which will be reported in detail else- 
where Jl^l) reveals that eigenstates which belong to the 
classical chaotic region are ergodic, in the sense that 
they fill a range of unperturbed energies, in a way which 



corresponds to the classical microcanonical distribution. 
Here we concentrate on the distribution of the occupation 
numbers of single-particle states. 

The distribution n s of occupation numbers of single 
particle levels s can be directly obtained from cigenfunc- 
tions. Given an eigenfunction ip n (E m ) one can write: 

n s (E m ) = ^ \ME m )\ 2 (n\h s \n) (3) 

n 

where n s is the occupation number operator. The term 
(n|n s |n) equals 0, 1, 2 depending on how many particles 
are located on the specific single-particle level s. In Fig.|l| 
some examples for the occupation numbers distribution 
(histograms on the left column) are given, together with 
the corresponding eigenstates (right column) . 
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FIG. 1. Eigenfunctions (right column) and the corre- 
sponding occupation number distributions (left column) for 
the case L — 3.5, / = 19. (ai) and (61) are for the 
ground state, (02) and (62) are for the 10* state (classically 
quasi-integrable), (03) and (63) are for the 49 th state (with a 
chaotic phase space and e m = E m /E ma x = —0.464) and the 
full line is the best fit to the BE distribution with j3fu = 0.305; 
(ai) and (64) are for the 54* state (classically chaotic) with 
e m = -0.449, P fit = 0.310 
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One can see a strong difference between the distribu- 
tions obtained from "regular" eigenstates, belonging to 
the stable region (61,62), and those obtained from er- 
godic ones, belonging to the chaotic region (63,64). In 
the stable region close eigenstates in energy E m yield 
completely different distributions n s (see Figjl] (a\, 0.2)), 
which means that there is no equilibrium in the statistical 
sense JTo| ] . 

In contrast, in the classically chaotic region the form of 
the distribution n s is statistically stable with respect to 
the choice of a specific chaotic eigenstate (compare Fig.[l] 
a 3 , 04). 

In the standard thermodynamical treatment, temper- 
ature can be defined in a number of different ways, which 
are known to be equivalent in the thermodynamical limit. 
It is interesting to compute temperature for our system, 
following these different definitions, although our system 
is far from this limit. First of all we use the canonical 
expression: 

y E -i3 a E m 
We. - y e ^E m (4) 

where E m are the exact eigenenergies of the interacting 
system. The above relation between energy and tempera- 
ture allows for standard thermodynamical description of 
our system. The solution (3 C (E) of equation (^) is shown 
in Fig.|| as a full curve. 

For a system of N non-interacting particles with to- 
tal energy E, the canonical distribution is well-known to 
result, when N is large, in the BE distribution for the 
occupation numbers: nf E = [e^ ri ( s+i )~^) — (s = 

—I, . . . , I) where \i is the chemical potential and j3 is the 
inverse temperature. This expression is derived for an 
ideal gas (many non-interacting particles) in contact with 
a thermostat; in contrast, our system is isolated, with two 
interacting particles only. Nevertheless, recent analyti- 
cal and numerical studies for random two-body interac- 
tion p|-|Io|] , suggest that conventional quantum statistics 
can appear even in isolated systems with relatively few 
particles, provided a proper renormalization of energy is 
taken. 

As a simple comparison, we may consider the BE dis- 
tribution as one-parameter fitting expression, taking into 
account the constraint set by the finite number of parti- 
cles (J2 S n f E = 2)- This allows to find the corresponding- 
inverse temperature (5 fa, and from that the BE distribu- 
tion presented in the left column of Fig.[l] by full curves. 

A different way of comparing numerical data with the 
BE distribution is to solve the following equations in the 
unknowns /3, fi: 

J2n? E = 2, J2 hsn s E =£ (5) 

s — — l s— — l 

with £ = hsn s , computed from numerical values of 
n s (note that £ is different from the exact eigenenergy be- 



cause our model is strongly non-perturbative p2[ ). Do- 
ing so, and averaging over a number of chaotic eigenstates 
with close values inside small energy windows, we have 
found inverse temperatures @be quite close to the previ- 
ously obtained [3fu (compare crosses with full circles in 
Fig.^). The agreement between the numerical values of 
(3fu and /3be supports, on one hand, the significance of 
the fitting procedure with a BE distribution and, on the 
other hand, the validity of the BE distribution for iso- 
lated systems with few interacting particles via a proper 
renormalization of the energy £, see [p|,[To| . 

We have then compared f3fu(E) with f3 c (E), where 
E is the exact eigenenergy, and have found them to be 
significantly different. This is of course hardly surpris- 
ing: even in the presence of ergodicity, with so few par- 
ticles, one cannot expect coincidence of microcanonical 
and canonical averages, which still strongly resent of the 
particular choice of the interaction. In classical terms, 
the distribution of single-particle energy still depends on 
the particular shape of the energy surface, which is in 
turn determined by the particular interaction chosen. 

A remarkably different result was obtained on "ran- 
domizing" our model, by replacing non zero off-diagonal 
matrix elements in the Hamiltonian by random variables 
with the same mean and variance as in the exact dynam- 
ical model (|l|). Whereas (3 C is left practically unchanged 
by this replacement, (3fu — (3be considerably changes, 
and it comes quite close to f3 c . In our understanding, the 
reason of this striking result is that the random model 
(which has no smooth classical limit) uses a much more 
"generic" form of the perturbation than the dynamical 
one; in addition, the dependence on the specific interac- 
tion is further weakened by the (matrix) ensemble averag- 
ing. As a result, our data indicate that the conventional 
canonical distribution of occupation numbers may appear 
even in an isolated system with a quite small number of 
particles with random interactions, where the dynamical 
correlations, which prevent a similar results in the clas- 
sical model, are negligible. 

Other different definitions of temperature have been 
widely discussed in application to complex nuclei [fl3| in 
the context of the onset of thermalization. It was in par- 
ticular shown that these definitions of temperature give 
the same result in realistic shell models of nuclei. This 
fact was connected with the onset of "true" thermaliza- 
tion. However, as was shown in fllPf , when the number of 
interacting particles is small, one can get different values 
of temperature depending on the particular definition, 
even in the equilibrium region, where a statistical de- 
scription is legitimate. In this connection we have also 
compared the temperatures found from the BE distribu- 
tion with the standard thermodynamical inverse temper- 
ature defined by, [3t = d Inp/dE where p is the density 
of states of the total Hamiltonian. Both for the dynam- 
ical and the random model, a Gaussian fit provides an 



3 



excellent approximation to actual data for the density 
of states 0. We have used this very fit in computing 
fir, obtaining practically the same result in both cases. 
The function (3t(E) is different from /3 C (E), but it comes 
closer and closer to it near the center of the spectrum 
(where both temperatures are infinite). For the dynam- 
ical model, (3t(E) turns out to be completely different 
from both (3 C and (3f it ~ (3be (see dotted line in Fig.||). 
This means that for small number of particles the above 
definition of (3t is irrelevant to the distribution of occu- 
pation numbers JTc| ] . For the random model, /3t is not 
far from /3fu(E) at high temperatures (where (3t and f3 c 
tend to coincide), but deviates from it at smaller tem- 
peratures. We do not push our comparison of 0fit(E) 
and p c to smaller temperatures than shown in Fig.gj, be- 
cause near the edges of the spectrum, where the density 
of states is small, the eigenstates are not fully chaotic 
any more (note that in the dynamical model the motion 
becomes increasingly regular there). 



generic fluctuations of n s for small changes of the energy 
thus invalidating any statistical approach. For chaotic 
eigenstates, the distribution of occupation numbers can 
be approximately described by the Bose-Einstein distri- 
bution, although the system is isolated and consists of 
two particles only. In this case a strong enough interac- 
tion plays the role of a heat bath, thus leading to thermal- 
ization. In spite of this, the minimal number of particles 
prevents the canonical distribution from describing our 
dynamical system even if a surprising agreement with the 
canonical distribution is recovered in the corresponding 
random model. 
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the University of Milan at Como for the hospitality dur- 
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FIG. 2. Dependence of the different definitions of inverse 
temperature /3 on the rescaled energy e = E/E max ; full curve: 



j3 c , dotted curve : f3r, full circles : j3fi 
circles are for the random model. 



crosses : (5b e- Open 



In conclusion, we have studied a dynamical model of 
two-interacting Boson particles in a finite dimensional 
Hilbert space. We have shown that, in the shell model 
representation, the structure of exact eigenstates can be 
directly related to the onset of equilibrium for the occu- 
pation numbers n s of single-particle states. Specifically, 
for the eigenstates corresponding to classical chaotic mo- 
tion, an equilibrium distribution for n s occurs which al- 
lows for a statistical description of the model. In con- 
trast, "regular" eigenstates results in extremely non- 
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